Nature of Intermediate States between Superfiuid and Mott insulator for Interacting 
Bosons in One-dimension with a Harmonic Trapping Potential 
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Successive quantum transitions in an intermediate regime are shown to exist between the super- 
fluid and Mott insulating states for interacting bosonic atoms in one dimension with a trapping 
potential. These transitions, which are caused by the interplay of the trapping potential with the 
competition between the kinetic energy and the interaction, reveal novel many-body effects as re- 
flected by low-lying excitation behavior, unconventional long-range correlations and an even-odd 
alternating squeezing process of superfiuid bosons into the Mott insulating state. These features, 
most likely being generic for all dimensions when a trapping potential is involved, are relevant for 
both experimental observations and physical interpretation of the Mott insulator transition. 
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Rapid developments in ultracold atom experiments 
on interacting bosons p| Q, fermionsj^, 0], and their 
mixtures 0, @] in recent years have greatly stimulated 
a full exploration of various fundamental properties of 
strongly correlated systems. One of the most important 
progress has been made with the experimental obser- 
vation of the Mott insulator transition from the super- 
fluid phase of interacting bosonic atoms in optical lat- 
tices with a harmonic trapping potential 0. The bosonic 
Mott insulator mimics the conventional Mott insula- 
tor with electron correlations, which has been challeng- 
ing condensed matter physicists for several decades. It 
was shown theoretically that bosons are uniformly dis- 
tributed and the Mott insulator transition occurs only 
at sufficiently strong interactions 0, H, @, HH- In 



experiments [2|, |12| . Il3l |. when a trapping potential is in- 



troduced to confine ultracold atoms, a shell structure for 
the local density shows up with a phase separation of 
Mott insulating and superfiuid atoms [3]. 

Experimentally, the Mott insulating state is exhibited 
by the vanishing of the interference pattern and the ap- 
pearance of resonances in the excitation spectrum. The 
Mott transition occurs when the visibility of interfer- 
ence pattern and its derivatives changes sharply with re- 
spect to the lattice depth. The larger quantum fluctua- 
tions in one dimension (ID) than those in 3D are shown 
from the Bragg spectroscopy [IBj]. Most recently, the mi- 
crowave spectroscopy with atomic clock shifts which pre- 
cisely measures the fluctuation of the atomic number 
at each site allows to identify the superfiuid layers be- 
tween the Mott shells 14| , while the formation of a dis- 
tinct shell structure and its incompressibility are shown 
for the in-trap density distribution by spatial selective 
microwave transitions and spin flip collisions. Numeri- 
cal simulations, while verifying the shell structure of the 
local density lid] , show additionally a kink structure of 
the visibility [l7|. All these findings raise an interesting 



issue whether and how the critical behavior of the Mott 
transition of the homogenous case is affected by the trap- 
ping potential. It seems that the single critical point for 
the homogenous case is expanded into an intermediate 
regime owing to the trapping potential. In this regime, 
we will show that those bosons in the central superfiuid 
region are compressed one by one into the Mott insulating 
plateau with increasing the interaction such that a series 
of successive quantum transitions take place due to the 
interplay of the trapping potential with the competition 
between the kinetic energy and interaction. 

In this study, we focus on exploring the nature of the 
intermediate states between the perfect superfiuid and 
Mott insulating states. For this purpose, we investigate 
the one-dimensional Bose-Hubbard model with a har- 
monic trapping potential, which can be realized properly 
by interacting bosons trapped in an optical lattice with 
one dominant recoil energy in experiments. The Hamil- 
tonian can be generally written as0, Q 
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H = -i^(oto i+ i + at +1 aiJ + ^-^ni(n*-l) 

i i 

+V T ^n l [ i -(L + l)/2] 2 (1) 

i 

where the lattice spacing has been set to one, t a hop- 
ping integral set as the energy unit, Vt the strength of the 
trapping potential with a center at x c = (L + l)/2, and 
U the on-site interaction. Here we employ density ma- 
trix renormalization group (DMRG) methodd M, [H 
which has proved very successful for studying several 
kinds of physical properties in quasi-lD strongly corre- 
lated systems. For the present system, we modify the 
finite-size algorithm in order to obtain a given accuracy 
for low-lying energy states and recover the spatial inver- 
sion symmetry with the trapping potential under open 
boundary conditions. The sweeping is conducted from 
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the middle to the two ends gradually rather than imme- 
diately and the convergence is reached for each sweeping 
length I G [4, L] and three lowest eigenstates are targeted 
simultaneously. Accuracies were examined for different 
sizes L, bare states per site n s , the number of states kept 
in the blocks m, and the total number of bosons TV and 
U. We found that n s = 8 and m = 100 are sufficiently 
large for U G [4,14], whereas n s = 32 and m = 400 
are necessary for U < 4 so that systematic errors due 
to truncations could be smaller than the symbol size (see 
below). While systematical studies were made for several 
properties with N G [20,70], novel effects in a crossover 
regime explicitly show up when TV > 50 at Vt — 0.01. In 
the following, our discussions focus on results for N = 60 
with Vt = 0.01, for which L = 80 is long enough to 
remove the boundary effects. 

Let us first discuss the low energy properties. For the 
homogenous case, it is well-known that the Mott tran- 
sition from the superfluid phase takes place at U c = 
3.61 ± 0.130 and a gap opens monotonically with re- 
spect to U in the Mott insulator phase. Fig. [1] shows 
gaps Ai.2 = -Ei. 2 — Eq for the two lowest excitations at 
Vt = 0.01 with So,i,2 being the three lowest energies. 
The phase diagram of the ground state now consists of 
three parts instead of two for Vr = 0. The first part 
is a weak coupling regime for U < U w including U = 
with U w « U c . In this regime (see inset), Ai slowly 
decreases until the emergence of oscillation at U = U w , 
while A2 behaves similarly except for a sharper change 
around U = 0. Since one has N = 60 harmonic oscil- 
lators with the center at x c for U = 0, all low-energy 
properties can also be obtained from numerical exact di- 
agonalization and perturbation calculations for small-C/. 
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FIG. 1: (color online). Energy gaps Ai,2 for the two lowest 
excitations relative to the ground state. A large-C/ regime for 
U > 12.12. Inset: for U G [0,4] at the same scale. 

In the second part, an intermediate regime with U w < 
U < U s where U s = 12.12, Ai^ vary smoothly at first 
and then behave as sawtooth with the emergence of the 
first subpeak. This behavior might be a consequence 
of the competition between the confining potential and 



the interaction. The peaks and dips for sawtooth re- 
sult from the level crossing. With further increasing of 
U, the ground state becomes degenerate for those small 
intervals of subpeaks, while E 2 can be equal to cither 
E\ or higher excitation energies in the intervals of main 
peaks. The dips indicate successive quantum transitions 
in the ground state. One has two kinds of excitations 
for two sides of each sawtooth (for both main and sub 
peaks). The left excitation is given by creating one dou- 
ble occupancy from the ground state, while the right one 
corresponds to the destruction of a double occupancy. 
Furthermore, the large-C/ regime, the third part with 
U > U s , is characterized by a unique ground state and 
is related to the fixed point of the large-C/ limit. At this 
fixed point, Eq is simply determined from the trapping 
potential and the local density rii, since no double occu- 
pancy is allowed and rii is simply divided into an inner 
part of n j = 1 , an outer part of rii = and an intervenient 
part with rii < 1. One therefore can consider this regime 
corresponding to the Mott insulator phase of Vt = 0. In 
the large-C/ regime, two lowest excitations are degener- 
ate. Due to the level crossing at U = 12.75, however, 
one has Ai )2 = U - 12.12 + 0(l/U 2 ) at first and then a 
Ai.2(= : 0.597) almost independent of U. The excitations 
are generated by creating a double occupancy from the 
ground state for 12.75 > U > U s and modifying the inter- 
venient part of the configuration of the ground state for 
U > 12.75. In addition, we have checked the low-energy 
behavior for odd and sufficiently large N, e.g. N = 59 
and found essentially the same features for U < U s but 
a two- fold degenerate ground state for U > U s . The de- 
generacy can be understood from the even-odd effect of 
the number of bosons in the large-C/ limit. 
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FIG. 2: (color online). The profile of m versus U and i. A 
staircase structure develops with increasing U and the upper 
part with Hi > 1 is eventually suppressed into a Mott plateau 
with rii = 1, i.e. a green interface. 

To interpret the above results, we show in Fig. [2] for 
U G [4, 14] how the profile of rii changes from a nearly 
global bell-shape at U = 4 into a complete plateau (green 
interface) that characterizes a Mott insulator with m = 1 



3 



for both i G (18,63) and the large-C/ regime. While the 
upper part (with red and yellow colors) of the bell-shape 
with n, > 1 gradually squeezes, a Mott plateau (green) 
grows from two flanks of horizontally at the height 
rii = 1 in the intermediate-C/ regime. Wc found very 
surprisingly that the squeezing is discontinuous for suffi- 
ciently large U with accompanying staircase structure for 
rii. Each staircase involves regularly an oscillation with 
respect to i but has different numbers of peaks. One can 
clearly see that the last three staircases have one, two 
and three peaks, respectively. We note that when U is 
increased, rii is first reluctant to squeeze by forming a 
bundle and then suddenly changes into its next bundle 
with one less peak at a critical value of U; The area of 
the upper part for each U is found precisely equal to the 
number of bosons involved in the upper part and the area 
between two-successive bundles is one within 5 digits ac- 
curacy, ft turns out that each peak essentially reflects the 
contribution of one extra boson to the central superfluid 
region. Therefore we can conclude that the superfluid 
bosons are squeezed one by one into the Mott plateau in 
the two flanks of nj when U is increased. 

1.0r 




FIG. 3: (color online). Bundles for single nf and double nf 
occupancies with an increment 0.1 for successive curves. 

To further elucidate the intrinsic physics for the 
squeezing process, we analyzed the contributions from 
different number occupancies to m for each site i, i.e. 
rii = J2 a n f where a = s,d,t, etc. nf, nf and n\ de- 
note single and double, triple occupancies of bosons at 
site i, respectively, nf can be obtained directly from the 
eigenvalues of the reduced density matrix pi at i. For 
U G [4, 14] , we found that nf and nf are larger than 
the rest by two orders of magnitude at least. Figure. [3] 
shows nf and nf only with U G [7.6, 14] for clearness. 
One can see that an oscillating structure similar to n, 
occurs for both nf and nf, but valleys in nf correspond 
to peaks in nf one by one. This implies the correlation 
between the single and double occupancies. Moreover, 



since the quantum fluctuation is not fully eliminated un- 
til the size of the Mott plateau £m is sufficiently large, 
one has nf = 1 — < 1 and n f = $(£m) f° r the large- 

U regime and i G (18, 63), although all superfluid bosons 
in this regime are squeezed such that ni = 1 precisely. 
For U < 12.12 and the curves with the same number of 
peaks, we found that the area between nf and <5(£jvf) hi 
Fig. [3](a) is twice of the corresponding area in Fig. [3jb) 
between 1 — 6(£m) and nf such that the area between 
ni and the Mott plateau precisely equals to the num- 
ber of superfluid bosons. This reveals that the superfluid 
bosons have phase coherence in a "quantized" way one 
by one when U is increased. 




FIG. 4: (color online) Successive transitions for n d , V, S ma x, 
and FWHM induced by increasing U. Red and green stair- 
cases correspond to colored curves in Figs. [3{a) and (b). 

Figure 5] shows, for the intermediate and large-?/ 
regimes, the average double occupancy n d . the visi- 
bility V, and the maximum value S max and the full- 
wide half-maximum(FWHM) for the momentum distri- 
bution, which were used to study the Mott insulator 
transition [l7j. In the largc-C/ regime, n d and S m ax 
change much slower than V and FWHM. When U is 
sufficiently large in the intermediate regime, a series of 
staircases develop for all four quantities as those shown 
similarly for n, in Fig. [5J One can easily see that long 
and short staircases show up alternatively and corre- 
spond to even and odd number of superfluid bosons in 
the central region of n*, respectively. Comparing this to 
Fig. [TJ one can also see that the shorter staircases cor- 
respond to those degenerate ground states with smaller 
intervals. The longer red and shorter green staircases in 
Fig. [1] correspond to the red and green curves in Fig. 
[3] as well as the last main and sub-peaks in Fig. [TJ In 
fact, once the Mott plateau of rij emerges at n^ = 1 in 
the squeezing process, the number of bosons involved in 
the intervenient as well as inner parts is even and is no 
longer changed. It turns out that squeezing one super- 
fluid boson enlarges the size of the Mott plateau by one. 
This implies that the degeneracy of the ground state, i.e. 
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the alternating nature, depends crucially on whether the 
number of supcrfluid bosons is even or odd. Moreover, 
when the ground state is non-degenerate, it would evoke 
larger quantum fluctuation to squeeze a superfluid bo- 
son into the Mott plateau since a two-fold degenerate 
ground state is subsequently generated. Therefore, the 
alternating nature reflects different quantum fluctuations 
in the squeezing process and a quantum transition in the 
ground state occurs at the edge of each staircase. 
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FIG. 5: (color online). One-body correlation functions p i , for 
various j with i 6 [1,80]. Black, red, green and blue curves 
are for j — i, 15, 22, and 39, respectively, (a)-(c) and (e) for 
the even parity ground state, but (d) with the odd-parity. 

Finally, we turn to the one-body correlation functions 
p^j = (±|ajaj|±) for different U in the intermediate 
regime as shown in Fig. [5] |±) stand for ground states 
with an even or odd parity, respectively There is always 
an even-parity ground state with the present trapping 
potential. The odd parity ground state occurs only when 
the ground state is degenerate for those small intervals 
of U with sub-peaks as seen from Fig. [1] corresponding 
to the shorter staircases in Fig. Q] One can see the sig- 
nificantly different behavior of pf ■ for different U and 
j as well. As expected, pfj decay faster for j in the 
Mott plateau than in the other regions. The correlation 
range extends over a regime with a non-zero occupancy 
of bosons as seen from (a) and (b), but is substantially 
suppressed when the size of the Mott plateau is suffi- 
ciently large as shown in (c)-(e). Moreover, for those U at 
which Ai reaches its local minima (see FigQ]), pf 15 (red 
curves) for j = 15 show that the correlation is greatly en- 
hanced with i in the right side of the intervenient regime 
as shown in (a)-(d). The enhancement is the largest when 
the ground state is degenerate as seen in (c) and (d). 

In summary, we have shown that a series of successive 
quantum transitions exist in the intermediate regime be- 
tween the superfluid and Mott insulator phases for in- 



teracting bosons in one-dimensional optical lattice with 
a harmonic trapping potential. Those transitions reflect 
the quantized squeezing process of the interacting bosons 
into the Mott insulator phase, revealing novel many-body 
effects resulting from the interplay of the trapping poten- 
tial with the competition between the kinetic energy and 
the many-body interaction. These induced effects are rel- 
evant for the interpretation of the Mott transition to the 
trapped atomic systems in all dimensions, in principle, 
since the trapping potential introduces an independent 
additional physical parameter, calling for further experi- 
mental explorations. 
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